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The isospin and strangeness dimensions of the Equation of State are explored. RIA and 
the SIS200 accelerator at GSI will allow to explore these regions in compressed baryonic 
matter. 132 Sn+ 132 Sn and 100 Sn+ 100 Sn collisions as well as the excitation functions of 
K/tt, A/ 7i and the centrality dependence of charmonium suppression from the UrQMD 
and HSD transport models are presented and compared to data. Unambiguous proof 
for the creation of a 'novel phase of matter' from strangeness and charm yields is not in 
sight. 

"o . 

The properties of neutron-rich matter at high densities and temperatures are largely 
unknown. Collisions of heavy rare isotopes provide an excellent tool to explore this un- 
known region in the phase diagram of nuclear matter. The planned Rare Isotope Acceler- 
ator (RIA) (up to E Lah = 500 AMeV) and the SIS200 at GSI (up to £ Lab = 30 AGeV) are 
the upcoming experimental facilities that will access new energy density regions for these 
collisions allowing to investigate the Equation of State (EOS) of compressed hadronic 
matter in regions (2-10 p neu tron) or even in directions (n neu tron = n \) which have not been 
experimentally accessible before. 

Among the most interesting open questions is the isospin and strangeness dependency 
of the EOS, especially at high baryochemical potentials. This has been the subject of 
intensive research, as it is not only important in the study of radioactive nuclei but 
bears consequences on various important astrophysical issues as well [T]. The EOS of 
isospin asymmetric matter can be written as e(p, 5) = e(p, 0) + E sym (p)5 2 + 0(<5 4 ), where 
5 = (p n — Pp)/(p n + P P ) is the isospin asymmetry, e(p, 0) is the energy density in isospin 
symmetric nuclear matter, and E sym (p) is a parametrisation for the isospin dependence 
of the nuclear symmetry energy [2]. Parametrisations for E sym (p) fall into two main 
categories jSj: (I) E sym (p) rises with p, (II) E sym (p) falls with increasing density p. The 
latter case has the interesting feature that symmetric matter may become unstable at 
high densities. 

The (tt~ /tt + ) ratio is particularly sensitive to the high density behaviour of E sym 
This can be easily understood qualitatively within the A resonance model for pion pro- 
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duction [1]: Here first-chance independent nucleon-nucleon collisions [5] yield a primordial 
(n- /n + ) ratio given by (5iV 2 + NZ)/(5Z 2 + NZ) « {N/Z) 2 . Thus, the charged pion ratio 
yields information on the isospin asymmetry (N/Z) in the dense participant matter, which 
depends on the high density behaviour of E sym [3]. The {ti~ /-n + ) ratio can therefore be 
used to probe the EOS of neutron-rich nuclear matter. 
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Figure 1. Time evolution of the (n/p) (top) and (n /ti + ) (bottom) ratios in central 
(6 = 3 fm) 132 Sn+ 132 Sn and 100 Sn+ 100 Sn collisions at 2 AGeV (left panels) and 11 AGeV 
(right panels) as calculated within the UrQMD model. 



For our investigation of this effect the UrQMD model vl.2 |6|7j is applied to 132 Sn+ 132 Sn 
and 100 Sn+ 100 Sn collisions at 2 AGeV and 11 AGeV beam energy. Fig. shows the 
time evolution of the (n/p) and (tt~ /ti + ) ratios. For 100 Sn+ 100 Sn, both the (n/p) and 
the (7r~/7r + ) ratios are roughly time independent and of order one, as expected for the 
symmetric system. For the asymmetric system 132 Sn+ 132 Sn, the (n/p) ratio drops during 
the compression phase from its initial value of 1.64 and then saturates. At 2 AGeV, 
this saturation occurs after 15 fm/c at a value of (n/p) fa 1.4, whereas at 11 AGeV the 
saturation occurs earlier (10 fm/c) and to a lower final value of (n/p) « 1.2. 

The (ti~ /-k + ) ratio rises slightly during the early phase of the reaction and saturates 
rather later than (n/p) (~ 25-30 fm/c at 2 AGeV and ~ 15-20 fm/c at 11 AGeV). The 
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saturation values of ~ 1.7 (2 AGeV) and ~ 1.2 (11 AGeV) fall far below the first-chance 
A resonance model prediction (N/Z) 2 « 2.7. This is evident, as it is not the initial but 
rather the local and dynamically changing (n/p) ratio that determines the (tt~ /tt + ) ratio. 
In addition, pion reabsorptions and rescatterings reduce the sensitivity of (n~ to 
(n/p). At higher energies, the increasing number of sequential nucleon-nucleon collisions 
drives the (ir~ /tt + ) ratio towards 1 via this mechanism. 

For isospin-asymmetric matter, a distillation effect has recently been proposed [S], which 
is powered by chemical instabilities in the liquid-gas phase transition driving the system 
towards isospin symmetry. This effect is similar to the well-known strangeness distillation 
mechanism in nuclear reactions, which has been proposed as a quark-gluon-plasma (QGP) 
signal 9J. The search for the QGP and its properties has been a driving force behind 
years of relativistic heavy ion physics. 

The fact that the A/tt ratio and the K + /tt + ratio are found experimentally to exhibit 
a maximum at -E , Lab=10-30 AGeV and drop to half that value at 160 AGeV has raised 
speculations about the appearance of the 'new phase transition'. The NA49 Collaboration 
has recently started an energy scan at the SPS. First results are available now at 40 and 
80 AGeV |10|11I12] and support these findings. Further studies are underway at 30 and 
20 AGeV H3|. 
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Figure 2. Excitation function of tc + and tt~ yields from central Au+Au (AGS) or Pb+Pb 
(SPS) collisions in comparison to experimental data |10|26j . Solid lines with open triangles 
show UrQMD results, solid lines with open squares and dashed lines stem from the HSD 
approach with and without potentials, respectively [2*3*] . 

Let us discuss these maxima in the excitation function on the basis of 'conventional' 
hadron transport theories - do they exhibit analogous maxima? In that case, the inter- 
pretation of the data as a QGP signal becomes doubtful. The various hadron spectra are 
conventionally calculated with non-equilibrium kinetic transport theory (cf. |14|15|16|17t 
I18|19j ). However, the calculated kaon to pion ratio from central nucleus-nucleus collisions 
turns out to vary by factors as large as 2 if different transport approaches are applied 
|19|20|21|22j . Thus a unique interpretation of the data is questionable so far. 
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Here we study the yields and ratios of protons, kaons, antikaons and hyperons from 
Au+Au (or Pb+Pb) collisions in the energy range from 1 AGeV to 160 AGeV [231 • The 
aim of this study is twofold: first, to find out the systematic differences between two 
currently used transport approaches (UrQMD |6l7j and HSD |14|24| ). and second, to 
compare to related experimental data |10|ll|12j . We provide predictions for experimental 
studies in the near future [JJJ], which are also of relevance for the new GSI-proposal [25J. 

The two different transport models employed, i.e. the UrQMD vl.3 (here an updated 
version has been applied) and HSD approach have been used to describe nucleus-nucleus 
collisions from SIS to SPS energies for several years. Though different in the numeri- 
cal realisation, both models are based on the same concepts: strings, quarks, diquarks 
(q, q, qq, qq) accompanied by hadronic degrees of freedom. It is important to stress that 
both approaches do not include any explicit phase transition to a quark-gluon plasma 
(QGP). A common failure of both models in comparison to experimental data may - 
model independently - indicate the appearance of a novel state of strongly interacting 
matter. 
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Figure 3. Excitation function of K + and K~ yields from central Au+Au (AGS) or Pb+Pb 
(SPS) collisions in comparison to experimental data [TO 26J. Solid lines with open triangles 
show UrQMD results, solid lines with open squares and dashed lines stem from the HSD 
approach with and without potentials, respectively E 



As a general overview on the 7T ± abundancies in central Au+Au and Pb+Pb collisions 
Fig. 121 shows the 7r + and n~ multiplicities at mid-rapidity as a function of the bombarding 
energy in comparison to the available data 10 2(i . Solid lines with open triangles show 
the results from the UrQMD calculations while the solid lines with open squares and 
dashed lines stem from the HSD approach with and without potentials, respectively. At 
lower AGS energies the UrQMD model gives slightly less pions then HSD (with/without 
potential), but both models overpredict the mid-rapidity data. Around 10 AGeV is the 
"crossing point" of both transport calculations and at SPS energies the tendency turns 
around: UrQMD gives more pions than HSD, so that HSD is now in a better agreement 
with the experimental data. 
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The K + and K~ multiplicities at mid-rapidity are depicted in Fig. El as a function of 
the beam energy in comparison to data from Refs. [10 26j: The K~ abundancies are well 
described by both transport models. The K + yield is slightly overestimated by UrQMD 
at AGS energies and underestimated at SPS energies, whereas HSD is in reasonable agree- 
ment with the data. Thus, Fig. El demonstrates that an underestimation of strangeness 
production is not the prevailing issue in comparison to the recent data from NA49 10 . 
Both transport models can roughly describe - within their systematic range of uncertain- 
ties - the X spectra and abundancies. 

Fig. 01 contrasts the K + /ir + and K~/tt~ ratios at mid-rapidity as a function of energy 
for central collisions of Au+Au (AGS) or Pb+Pb (SPS) with the data |10l2ti| . The 
excitation function of the K~ /7r~ ratio is roughly reproduced by both transport models, 
the maximum in the K + /n + ratio seen experimentally is not described quantitatively by 
either HSD or UrQMD. For the K + /tt + ratio both models give quite different results. 
HSD gives a monotonous increase of this ratio with bombarding energy (in qualitative 
disagreement with the data as pointed out in Refs. |20I21| ). whereas within UrQMD the 
ratio shows a maximum around 10 AGeV and then saturates to a slightly lower value 
towards RHIC energies, thus reproducing the behaviour of the data qualitatively. In view 
of Figs. [21 and 01 this failure is not primarily due to a mismatch of strangeness production, 
but more due to an insufficient description of the pion abundancies. 




(AGS) or Pb+Pb (SPS) collisions in comparison to experimental data |10|26j . Solid lines 
with open triangles show UrQMD results, solid lines with open squares and dashed lines 
stem from the HSD approach with and without potentials, respectively [2*3*| . 



This overestimation of pion abundancies can be attributed to a number of reasons: One 
problem of the transport approaches used here is that detailed balance is not implemented 
for n <-> m transitions with n, m > 3 |27j . Thus multi-particle collisions might change 
the dynamical picture accordingly and lead to 'shorter' chemical equilibration times [2H] ■ 
In fact, the importance of 3 <-> 2 transitions has been explored in an extended HSD 
transport approach [2E] for antiproton reproduction by meson fusion for A + A collisions 
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at the AGS and SPS. In order to achieve a more conclusive answer from transport studies 
multi-particle interactions deserve further investigation in future generations of transport 
codes. 

Another reason for the overestimation of the pion yield might be that the pions in 
both transport models are treated as 'free' particles, i.e. with their vacuum mass. Lat- 
tice QCD calculations and effective Lagrangian approaches like the Nambu-Jona-Lasinio 
(NJL) model however, indicate an increase in the pion mass with temperature and den- 
sity. So the overestimation of the pion yields in HSD and UrQMD might be a signature 
for a dynamically larger pion mass. Moreover, in-medium changes of the strange hadron 
properties, as known from experimental studies at SIS energies, may also show up in 
the compressed baryonic matter encountered in the AGS and SPS energy range. Thus, 
including all medium effects simultaneously in a consistent way might provide a more 
conclusive interpretation of the ratios in Figs. 0] and El However, such calculations re- 
quire a precise knowledge about the momentum and density dependence of the hadron 
self-energies which is not available so far. 




# exp. data 
—□—HSD, pot. 

HSD, cas. 

—A— UrQMD 



10 

E, VA [GeV] 



100 



0.08 



© 
II 

>> 



0.06 . 



0.04 



0.02 



0.00 




(A+2L ) / k 



# exp. data 
—□—HSD, pot. 

HSD, cas. 

—A— UrQMD 




10 

E lah /A [GeV] 



100 



Figure 5. Excitation function of the A + £° yields and (A + S°)/7r ratio from central 
Au+Au (AGS) or Pb+Pb (SPS) collisions in comparison to experimental data [10 26J. 
Solid lines with open triangles show UrQMD results, solid lines with open squares and 
dashed lines stem from the HSD approach with and without potentials, respectively [23J. 



Fig. El shows the excitation functions of A + S°-hyperons (left) at mid-rapidity as a 
function of the bombarding energy for central collisions of Au+Au (AGS) or Pb+Pb 
(SPS) in comparison to data jllll2|29|3U|31j . Here both models compare rather well 
with data. The (A + S°)/7r ratios (right) at mid-rapidity are underestimated slightly 
which should again be attributed to the pion excess in the transport models (as discussed 
above). Nevertheless, the experimentally observed maxima in the ratios are qualitatively 
reproduced by both models, indicating that with increasing bombarding energy s-quarks 
are more frequently produced within mesons (K, K ) rather than in associate production 
with baryons. A similar trend is also found in statistical models |32j . 

What can we conclude from the common failures of the models studied in comparison 
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Figure 6. The ratio of J/^> to Drell-Yan production as a function of forward energy Ezbc 
(left) and transverse energy (right) in Pb+Pb collisions at 158 AGeV. Data points are 
preliminary results from the year 2000 NA50 Pb+Pb run (Hill, the histograms show the 
1999 UrQMD results [H]. 

to related experimental data and does this failure provide evidence for a new state of 
matter in view of a QGP? In view of the large 'systematic uncertainties' in the discussed 
transport approaches we can quote an insufficient accuracy in the description of the pion 
degrees of freedom by both transport models. Thus, the deviations between theory and 
the data should not be seen as unambiguous signals for any kind of new physics. 

Another much-discussed possible QGP signal is charmonium production. Here, new 
data on the end Ezbc dependence of J/\I/ production in Pb+Pb collisions at 158 AGeV 
from the year 2000 NA50 run has recently become available - This data (one analysis), 
along with UrQMD results [33] (which are in line with HSD results [33]) is shown in Fig. El 
Note that the puzzling drop of the J/^> to Drell-Yan production at high Et which was 
present in the 'old', previously available NA50 data is no longer present in the current 
NA50 data. The new NA50 data fall very close to the published UrQMD calculations by 
Spieles et al. [3J]- To our opinion, both the £+ dependence and the Ezdc dependence do 
no longer justify to speak of an 'anomalous' J/^f suppression in high energy heavy ion 
collisions. 

We have shown that the exploration of the isospin and strangeness dimensions of the 
Equation of State opens fascinating new areas of research. Especially the upcoming 
experimental possibilities like the planned Rare Isotope Accelerator RIA and the SIS200 
machine at GSI will allow to explore novel phenomena in compressed baryonic matter in 
great detail. While conclusive unambiguous proof for the creation of a 'novel phase of 
matter', e. g. Quark Gluon Plasma, from strangeness and charm yields is still lacking, the 
hunt and quest for it has provided new insights and even opened novel fields of research. 
It has become a driving force behind modern high energy nuclear physics. 
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